We present a compositional analysis of 8 - 13ji/m spec- 
tra of 32 young stellar objects (YSOs). Our sample consists 
of 5 intermediate-mass stars and 27 low-mass stars. While the 
spectra and first scientific results have already been published 
by Przygodda et al. (2003) and Kessler-Silacci et al. (2004) 
we perform a more detailed analysis of the ~ 10 //m silicate 
feature. In our analysis we assume that this emission feature 
can be represented by a linear superposition of the wavelength- 
dependent opacity K a bs(/1) describing the optical properties of 
silicate grains with different chemical composition, structure, 
and grain size. The determination of an adequate fitting equa- 
tion is another goal of this study. Using a restricted number 
of fitting parameters we investigate which silicate species are 
necessary for the compositional fitting. Particles with radii of 
^^ ' 0. 1 fim- and 1 .5 //m consisting of amorphous olivine and pyrox- 
f^i , ene, forsterite, enstatite, and quartz have been considered. Only 
CN compact, homogeneous dust grains have been used in the pre- 
!— i , sented fitting procedures. In this context we show that accept- 
£jH' able fitting results can also be achieved if emission properties 
■^ , of porous silicate grains are considered instead. Although some 
■^i- ' previous studies give reasons for the similarity between the dust 
in circumstellar disks of T Tauri stars and Herbig Ae/Be stars, 
a quantitative comparison has been missing, so far. Therefore, 
we conclude with a discussion of the results of a 1 fim spectro- 
scopic survey of van Boekel et al. (2005) who focus on Herbig 
f~ ■j , Ae/Be stars, the higher mass counterparts of T Tauri stars and 
■^J" draw comparisons to this and other studies. We find that the re- 
CD suits of our study of T Tauri systems partly agree with previous 
^^ ■ studies of Herbig Ae/Be stars. 

,-G Key words. Astrochemistry - Stars: pre-main sequence - 
r ~ ( Infrared: stars 
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Abstract. 



1. Introduction 

TTauri stars are known as precursors of low-mass main se- 
quence stars similar to our sun. In contrast to their more 
evolved counterparts one of their characteristics is excess ra- 
diation observed in the infrared and sub-millimeter range. 
High-resolution interferometric observations (e.g., Leinert et 
al. 2004; van Boekel et al. 2004; Millan-Gabet et al. 2001) com- 
bined with radiative transfer models (e.g., Cotera et al. 2001; 
Wolf et al. 2003c; Akeson et al. 2005) confirm what has been 
assumed for 30 years (Mendoza 1968; Adams et al. 1987): the 
excess arises from a circumstellar, passive disk heated by the 
central star. Last but not least images in the near-infrared up to 
the sub-millimeter wavelength regime have provided impres- 
sive illustrations of the circumstellar structure of young stellar 
objects (YSOs) (e.g., Padgett et al. 1999; Allen et al. 2002; 
Krist et al. 2005). 

Several different absorption and emission features have 
been detected towards YSOs, so far (e.g., Molster & Waters 
2003 and Table 1 in Voshchinnikov 2002). A broad, frequently 
found emission band is the silicate feature at ~ 10 fim which 
is attributed to the stretching mode of the Si-O bond in sili- 
cate minerals (e.g., Dorschneret al. 1988). Hanneret al. (1995) 
successfully reproduced this feature and the underlying contin- 
uum by a two component fit: XT a k CqA" + C\e{X)A m , where 
Co and C\ as well as m and n are fitting constants and e{X) 
is the normalised, smoothed silicate band of 8 Ori D (Elias 
1978) in the N band '. The first term represents the underlying 
continuum from the optically thick, inner region of the disk, 
while the second term stands for the emission feature which 
is assumed to have its origin in the optically thin surface layer 
(Natta et al. 2000). Instead of using the normalised emission 
feature of an arbitrarily selected source, Bouwman et al. (2001) 
performed compositional fits which are based on the emission 
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profiles of different silicate compounds, particle shapes and 
radii (0.1/im and 2.0//m for olivine, 0.1/im for forsterite, en- 
statite and quartz). This approach was motivated by previous 
studies of the optical parameters of different silicate minerals 
(e.g., Dorschner et al. 1995). The resulting profile is a linear 
combination of the opacities Kabs^W) f° r the different silicate 
species i. 

In respect of the observed shift in peak position of the 
silicate band, Bouwman et al. (2001) concluded a change in 
average grain size and dust crystallisation in Herbig Ae/Be 
(HAeBe) stars indicating dust grain evolution. In addition, this 
result was supported by further comparisons with mid-infrared 
spectra of the interstellar medium as the parent material for 
YSOs, and the mid-infrared spectra of the comet Hale Bopp 
as a representative of the evolved planetary systems. While the 
mid-infrared flux of interstellar medium shows a maximum at a 
wavelength of ~ 9.8^m, the silicate feature of comet Hale Bopp 
is flatter with a maximum at ~ 11.3/im. Corresponding results 
were found in laboratory for silicate grains of increasing size 
and crystallinity. 

Because of instrumental sensitivity constraints most studies 
have focused on HAeBe stars rather than on the fainter T Tauri 
systems up to the present. Apart from the similarity of general 
properties like geometry and structure, evolutionary processes 
like dust coagulation and crystallization are assumed to be ac- 
celerated or at least enabled in HAeBe stars because of the 
higher energy input from the star. There are only a few authors 
who analysed the 10//m-feature of classical TTauri stars, so 
far. Natta et al. (2000) presented low-resolution spectra of nine 
classical TTauri stars associated with the Chameleon I dark 
cloud. Honda et al. (2003) found a larger amount of crystalline 
silicate like forsterite, enstatite and quartz in the primary of the 
well-studied binary Hen 3-600. Furthermore, Przygodda et al. 
(2003, hereafter PR03) and Kessler-Silacci et al. (2005, here- 
after KS05) selected a sample of 16 TTauri and 5 HAeBe stars. 
They established a linear correlation between feature strength 
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and feature shape which were explained by grain growth. The 
similar correlation has been ascertained for HAeBe stars (van 
Boekel et al. 2003, 2005). 

In this paper we complete the work presented in PR03 
and KS05 by providing a more detailed analysis of the sili- 
cate emission feature. While the objects are introduced in Sect. 
2, in Sect. 3 the fitting function is scrutinised in detail in re- 
spect of the temperature (Sect. 3.1) and the emission properties 
(Sect. 3.2) of the dust grains which contribute to the silicate 
feature. After studying our fit results (Sect. 4) we investigate 
the dependence of grain growth and crystallization on stellar 
parameters and look for further correlations. Finally, we raise 
the question, if the silicate feature could be reproduced by dust 
species that have not been considered, so far (Sect. 5). In this 
context, porous dust grains are taken into account. A discussion 
of future prospects conclude this study in Sect. 6. 



In fact, our first approach was the determination of a poten- 
tial grain size distribution of dust in a circumstellar disk, only 
considering the 10//m silicate feature. This approach failed be- 
cause unique fitting models could not be obtained. In this con- 
text, we have to mention that Mathis et al. (1977) found a par- 
ticle size distribution for graphite and several silicates in in- 
terstellar dust that can roughly be reproduced by a power law. 
It is still an open issue if a similar particle size distribution is 
existent in the environment of YSOs, too (e.g., Tanaka et al. 
1996). 

Nevertheless, it has been found, that a simple function with 
only a small number of fitting parameters finally guarantees a 
successful, robust fit. Therefore, the balance between physical 
justification and fitting robustness is our guiding idea. 



2. The sample 

In Table 1, our sample of 27 TTauri is presented. We addi- 
tionally consider 5 HAeBe stars in order to show differences 
and similarities between low-mass and intermediate-mass sys- 
tems. Furthermore, we are able to analyse differences between 
our results and previous studies. The targets were observed in 
the N band with the Thermal Infrared Multi Mode Instrument 
2 (TIMMI2; Reimann et al. 1998, 2000) at ESO's observatory 
La Silla and with the Long Wavelength Spectrometer (LWS) at 
the W.M. Keck Observatory (e.g., Marler et al. 1995), respec- 
tively. Information about the observations and the subsequent 
data reduction are given in PR03 and KS05. 

Six of our targets are known binaries (e.g., VW Cha, Hen 3- 
600 A, AK Sco), but we have seperate spectra only for two bi- 
nary objects (AS 205 and S CrA). Therefore, we present the 
spectral analysis of (30+2) objects altogether. Selection ef- 
fects in terms of stellar age and stellar mass can not be found 
(see Table 1), although most of our objects are younger (age 
< 3Myr) TTauri stars. Assuming similar stellar masses, the 
putative more frequent observation of younger T Tauri stars can 
be explained by their stronger infrared brightness in contrast 
to older TTauri stars (> 6Myr) where the circumstellar disk, 
which is responsible for the infrared excess, has already van- 
ished (Haisch et al. 2001; Carpenter et al. 2005). The age is 
certainly the least constrained parameter in Table 1 . 

Considering the inclination of the circumstellar disks of the 
objects as a further, geometrical selection effect, the appearance 
of an emission feature favours smaller inclination angles 2 if the 
feature originates from the innermost regions of the optical thin 
surface of the disk. 



3. The fitting function 

The determination of an adequate fitting function of the 10 pm 
silicate feature and its physical justification is non-trivial. The 
difficulty results from the attempt to simulate the complex spec- 
tral energy distribution (SED) with the least number of para- 
meters. 



with respect of the face-on orientation 



Table 1. Our sample of 27 TTauri stars and 5 HAeBe stars (lower 5 panels) with the corresponding coordinates. 





# 


Object 


RA (J2000) 


DEC (J2000) 


7*[K] 


M+ [M Q ] 


c/[pc] 


Sp. Type 


L+ [L ] 


age [Myr] 


Other common names 


Instrument 




1 


GGTau 


04 32 30.3 


+ 17 3141 


38020) 


1.3( 3 ) 


140 


K7<2) 


1.6(3) 


0.7-1.7<3) 




TIMMI2 




2 


AATau 


04 34 55.6 


+24 28 54 


4060( 4 ) 


0.7(5) 


140(4) 


K7<5) 


1.0*5) 


2.4(4) 




LWS 




3 


LkCa 15 


04 39 17.8 


+22 21 05 


4350< 4 > 


l.O* 3 ) 


140<4) 


K5 V<4) 


0.7(3.4) 


2.0( 4 ) 




LWS 




4 


DQTau 


04 46 53 


+ 17 00 00 


4000( 21 ) 


0.4(22) 


140(20 


M0< 7 > 


0.6(7.22) 


- 




TIMMI2 




5 


DRTau 


04 47 05.5 


+ 16 58 42 


4060< 3 - 6 > 


1.1(6) 


140(D 


K4V 


1.0(3.7) 


2.5-3.8(3) 




TIMMI2 




6 


GMAur 


04 55 10.9 


+30 22 01 


4060( 4 ) 


1.2(3.6) 


140<4) 


K7(4) 


1.0*4) 


1.3-1.8(3) 




LWS 




7 


SUAur 


04 55 59.4 


+30 34 02 


5945( 8 ) 


1.7(8) 


152 


Gl<8) 


7.8<8) 


8.7(8) 


HD 282624 


TIMMI2 




8 


GWOri 


05 29 08.4 


+ 1152 13 


6030( 8 ) 


3.7(8) 


305 


K3 W> 


61.8( 8 ) 


1.0(8) 


HD 244138 


TIMMI2 




9 


FUOri 


05 45 22.6 


+09 04 12 


5000( 23 ) 


1.0(23) 


- 


G3I 


- 


- 




TIMMI2 


B 


10 


BBW76 


07 50 35.5 


-33 06 24 


5500(24) 


- 


180(24) 


G0-G2 1(24) 


- 


- 




TIMMI2 


■■> 

eh 


11 


CRCha 


10 59 07.0 


-77 01 40 


4900 


1.2C) 


144 


K2 


3.3(9) 


1.0(8) 




TIMMI2 


| 


12 


TWHya 


110151.9 


-34 42 17 


4000< 10 > 


0.6< 10 > 


55< 10 > 


K8Ve 


0.2(1°) 


10(1°) 




TIMMI2 


<* 


13 


VWCha 


1108 01.8 


-77 42 29 


4350P) 


0.6(9) 


- 


K2/K5<9> 


2.9(9.26) 


1.0(9) 




TIMMI2 


c 
u 


14 


Glass I 


1108 15.4 


-77 33 53 


4600O 


0.9<9) 


- 


K4<u> 


1.6*9) 


1.0(9) 




TIMMI2 


t 


15 


VZCha 


1109 23.8 


-76 23 21 


4200(25) 


0.7*25) 


- 


K7<25) 


0.5(9.25) 


3.0(9) 




TIMMI2 


O 

'1 


16 


WWCha 


11 10 00.7 


-76 34 59 


4350(H) 


0.6<"> 


- 


K5 


2.7(H) 


0.3(H) 




TIMMI2 


g 


17 


Hen 3-600 A 


11 10 28.9 


-37 32 05 


3200( 5 ) 


0.3C2) 


50(4) 


M4<4) 


0.2(12) 


10(4) 




LWS 


o 
u 

r—. 


18 


IRAS 14050-4109 


14 08 10.3 


-41 23 53 


4405(4) 


- 


140(4) 


K5<4) 


0.6*5) 


10<5) 




LWS 


y: 

B 


19 


Sz82 


15 56 42.3 


-37 56 06 


3800d5) 


0.3(26) 


- 


M0<26) 


1.6(26) 


1.1(26) 


IMLup 


TIMMI2 


3. 

C 


20 


RULup 


15 56 42.3 


-37 49 16 


4000< 13 > 


0.8C3) 


140< 14 > 


G5 V 


0.5(13) 


0.04d5) 


HD 142560 


TIMMI2 


•3 


21 


AS 205 N 


16 1131.4 


-18 38 25 


4450< 16 ) 


1.5< 16 > 


160d6) 


K5<i6> 


7.1(16) 


O.K'6) 




TIMMI2 




22 


AS 205 S 






3450< 16 > 


0.3C6) 


160( 16 > 


M3<i6> 


2.2(16) 


<0.1<16) 




TIMMI2 




23 


HBC639 


16 26 23.4 


-24 21 02 


5250< 10 ) 


2.3(io) 


160( 10 > 


K3<i°> 


8.8(1°) 


1.5(1°) 


DoAr24 


TIMMI2 


o 

GO 


24 


Harol-16 


16 3133.5 


-36 53 19 


4365< 15 > 


- 


- 


K2-3d 7 > 


2.0d5) 


0.5(15) 


DoAr44 


TIMMI2 




25 


AKSco 


16 54 44.9 


-36 53 19 


6500< 18 > 


- 


152d8) 


F5V 


0.9d9) 


6.0d9) 


HD 152404 


TIMMI2 




26 


SCrAN 


19 0108.6 


-36 57 20 


4800< 16 > 


1.5d6) 


130d 6 ) 


K3<i6> 


2.3(16) 


3.0d6) 




TIMMI2 




27 


SCrAS 






3800< 16 > 


0.6C6) 


130d6) 


M0 


0.76d6) 


1.0(16) 




TIMMI2 




28 


MWC 480 


04 58 46.3 


+29 50 37 


8670( 4 ) 


- 


145<4) 


A3<4) 


22<4) 


4.6 


HD 31648 


LWS 




29 


HD 163296 


17 56 21.3 


+29 50 37 


9332(4) 


2.0(20) 


122<4) 


Al V<4) 


24.0<2°) 


4.0<4) 




LWS 




30 


HD 179218 


19 11 11.3 


15 47 16 


10471(4) 


2.9(20) 


240<4) 


B9<4) 


100<2°) 


0.1<4) 


MWC 614 


LWS 




31 


WWVul 


19 25 58.8 


+21 12 31 


8600(4) 


- 


550(4) 


A3/B9 V<4) 


- 


- 


HD 344361 


LWS 




32 


HD 184761 


19 34 59.0 


+27 13 31 


7500(4) 


- 


65<4) 


A8V<4) 


7.4(4) 


- 




LWS 



Note: Following symbols are used: stellar temperature (7"*), stellar mass (M*) derived from corresponding stellar evolutionary tracks, distance (<af), stellar luminosity (L+), and stellar age (mostly derived from their location in the HR diagram). 
References - (1) Thi et al. 2001; (2) Hartmann et al. 1998; (3) Greaves 2004; (4) KS05; (5) Metchev et al. 2004; (6) Kitamura et al. 2002; (7) Muzerolle et al. 2003; (8) Calvet et al. 2004; (9) Meeus et al. 2003; (10) Natta et al. 2004; (11) Natta et 
al. 2002; (12) Geoffray & Monin 2001; (13) Stempels & Piskunov 2003; (14) Brooks & Costa 2003; (15) Niirnberger et al. 1997; (16) Prato et al. 2003; (17) Chen etal. 1995; (18) Guntheretal. 2004; (19) Mansetet al. 2005; (20) van Boekel etal. 2005; 
(21) Mathieu et al. 1997; (22) Hersant et al. (2005); (23) Lachaume et al. 2003; (24) Schiitz et al. 2005; (25) Natta et al. 2000; (26) Gomez & Mardonez2003. 
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In order to compare our results with previous ones, we con- 
sider an ansatz which Hanner et al. (1995) and BO05 used: 



r v 



> . CiK* 



By(T)\C + / / ,A' a hs:, 

!=1 



(1) 



where Co and C, are the fitting parameters, k^-j represents the 
wavelength-dependent opacity for dust component i, and B Y (T) 
is the Planck function corresponding to the temperature T (see 
Eq. 2). In the Sects. 3.1 and 3.2 we will include more detailed 
discussions of Eq. (1). 

Two fitting procedures are considered. First, we examine 
a fitting procedure that is based on an iterative Levenberg- 
Marquardt least-squares minimization (Press et al. 1986). But, 
this numerical fitting method does not guarantee the best fit, i.e. 
the finding of the global minimum of the fit quality parameter 
X 2 (reduced chi-squared, see Sect. 4), particularly for data with 
lower signal-to-noise ratios. Repeated modifications of the ini- 
tial fitting parameters and the subsequent analysis of x 2 allows 
an effective search for the global minimum of x 2 - 

Second, we consider the non-negative, linear least-squares 
fit procedure used by Lawson & Hanson (1974) containing the 
analytical solution of the matrix equation CjK a bs;i = fy with 
£abs;i = ByKubsj- The matrix equation is solved for each integer 
temperature value between a minimum and maximum temper- 
ature. The temperature with the lowest corresponding^ 2 is the 
temperature T that we look for. We assume a silicate sublima- 
tion temperature of 1500 K (Duschl et al. 1996) and a negli- 
gible contribution of the silicate emission to the mid-infrared 
spectrum if the temperature of the corresponding silicate com- 
pound is lower than 50 K (see Sect. 3.1). Because of robustness 
and high rapidity the non-negative, linear least-squares fit pro- 
cedure is finally used for our study. 

3.1. The single-temperature black body approximation 

In Eq. 1 the Planck function has the following form (x = 
(Av)/(*T)): 



15 

B V (X) = -J 



iv exp(x) - 1 



(2) 



with the Planck constant h, the Boltzmann constant k, and the 
frequency v. This function is normalised in such a way that its 
integral L B v (T)dv is unity. 

In fact, dust grains with different temperatures contribute to 
the SED in the N band. But an extensive study of the silicate 
feature and the underlying continuum can only be done with 
radiative transfer simulations in a dust disk model. However, 
there have been several different efforts in order to determine 
the different silicate compounds with simple, semi-analytical 
tools. Most models are based on the idea that the silicate feature 
arises from single-temperature dust grains in an optically thin 
surface layer. Consequently, the emitted flux is given by the 
formula T v - t v B v (T) with the optical depth t v = £" =1 C,-/c a bs;i- 
Although this approach does not consider the dependence of 
dust temperature on grain size and composition, the mean dust 
temperature is assumed to be dominated by the absorption 
and emission characteristics of carbon grains assuming thermal 



equilibrium (Kriigel & Walmsley 1984). Przygodda (2004) ob- 
tained a further simplification of the model by the assumption 
that the temperature of most dust grains, which effectively con- 
tribute to the emission feature, is T ~ 300 K. As a maximum 
flux [B v (T)] m - dx at A - 10;um corresponds to a temperature of 
~ 290 K the Planck function B V (T » 300 K) is approximately 
constant in N band and the shape of the feature depends only 
on the emission properties of dust grains: T v ~ T v . In order 
not to preselect a certain dust temperature, the temperature T 
can also be introduced as another free fitting parameter as it 
was done in BO05 and as we do in this study. Also in this ap- 
proach the underlying continuum of the 10//m-feature is as- 
sumed to be dominated by dust grains of the same tempera- 
ture regime. We have to mention that this approach does not 
consider the possibility that large reservoirs of cool (<k300K) 
and/or hot (»300 K) dust grains provide the main contribution 
to the SED at 10 fim. For this reason a further analysis requires 
the modeling of a more representative, circumstellar disk. Such 
a representative model of a circumstellar disk is generated by 
means of a Monte Carlo code for radiative transfer simulations 
(MC3D: Wolf et al. 1999; Wolf 2003a,b). The density distribu- 
tion is calculated assuming hydrostatic equilibrium in vertical 
direction, i.e. perpendicular to the disk midplane (Schegerer et 
al., in preparation). 

For this study, stellar properties are chosen that were re- 
cently found by modeling the SED of the low-mass star RY Tau 
(Schegerer et al., in prep.; Akeson et al. 2005). The inner and 
outer radius of the disk amounts to 0.3 AU and 100AU, respec- 
tively. As absorbing, scattering and emitting material we use a 
mixture of carbon and silicate dust in a 1:2 ratio (Weingartner 
& Draine 2001; Draine & Lee 1984; Mafhis et al. 1977) 3 . This 
dust mixture has been used in a large variety of simulations 
of circumstellar dust configurations (e.g., Wolf et al. 2003d). 
The extinction, absorption and scattering cross sections are av- 
eraged over their abundance. Such an approach minimises the 
required computer power and memory as the temperature and 
density distribution has not to be determined and stored for 
each dust type (Wolf 2003c). Figure 1 shows spectra that are 
obtained from systems with identical geometry and identical 
stellar properties but in each case only dust grains having a 
temperature in the indicated temperature interval are taken into 
account. Although the 10 //m flux is maximised for a temper- 
ature T of about 375 K (4) , dust grains with a temperature be- 
tween 300 K and 400 K contribute only moderately (~ 20%) to 
the spectra compared with grains in adjacent temperature inter- 
vals. Therefore, the N band of the SED can not be assigned to 
a certain temperature interval or a certain area of the disk. In 
reality, all grains with different temperatures contribute to this 
feature. 

However, the temperature that we obtain in our single- 
temperature approximation is an average temperature of all 
dust grains weighted by their spectral contribution to the ob- 
served silicate feature. The average itself is assured by the 
monotonic gradient of the spectral contribution of dust grains 



3 See also http://www.astro.princeton.edu/"' draine. 

4 This value is obtained by a more detailed analysis with more nar- 



row temperature intervals. 
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temperature interval: 100. K - 200. K 



temperature interval: 200. K - 300. K 



temperature interval: 300. K - -100. K 









A [ M m] 

temperature interval: 800. K 



A h 
temperature interval: 900. K 






Fig. 1. The simulated SEDs of a T Tauri system assuming a distance of 140pc. The inner and outer radius of the disk amounts to 0.3 AU and 
100 AU, respectively (see Sect. 3.1). The spectra are obtained from systems with identical geometry and identical stellar properties but in 
each case only dust grains that have a temperature in the indicated temperature interval contribute to the SED. Each figure shows the SED for 
inclinations 0°, 15 ° and 45 ° with respect of face-on orientation. Each SED is calculated for 34 wavelengths (O). The integration of the flux 
between 5.0,um and 15.0yum yields 1.84, 3.9, 7.5, 5.1, 5.5, 4.3, 4.8, 2.8, 1.8xl0 14 Jy x Hz from the upper left to the lower right figure. The 
feature at ~ 20/iin is ascribed to silicate, too. 



with temperatures in the interval between 100K and 600 K. 
Figure 2 shows the fitting quality parameter x 1 ( see Eq. 3) of 
some fits versus the pre-specified temperature T for the Planck 
function (2). The existence of a global minimum of (x(T)) 2 jus- 
tifies the simplifying single-temperature-model of the underly- 
ing spectral continuum in N band. 

Future radiative transfer simulations will give the neces- 
sary assistance for the decision if our simplified model for 
silicate emission in circumstellar disk is physically reason- 
able. Furthermore, interferometric observation with the Mid- 
Infrared Interferometer (MIDI) could help to understand the 
potential, physical meaning of the determined, fitting tempera- 
ture T (Schegerer et al., in prep.). 

3.2. Emission profiles 

The emission profiles that we use for the fitting procedure 
are calculated for homogeneous, compact, spherical grains 
with Mie theory starting from the complex refractive indices 
n, for silicate material i (Dorschner et al. 1995; Servoin & 
Piriou 1973; Jiiger et al. 1998b; Spitzer & Kleinmann 1960). 
The result is the dimensionless absorption efficiency Q d bsj- 



This efficiency is used for calculations of the opacity 5 /e a b s; ; = 
2abs;if" 2 /(4/37rr 3 p), where r is the particle radius and p the 
material density. 

Fig. 3 shows the opacities for amorphous olivine 
(MgFeSiCU) and pyroxene (MgFe[SiC>3]2) while the lower 
panels present crystalline forsterite (Mg2SiC>4), enstatite 
(MgSiCh) and quartz (SiCh), each with radii 0. 1 fim, 1 .5 fim, 
and 2.5 //m. As Bouwman et al. (2001) pointed out, one 0.1/im- 
and one 2.0/im-sized dust component for each chemical com- 
pound seem to be sufficiently representative to model the 
10/im-silicate-feature. Here, we use dust grains with a radius 
of 0. Ifim and 1 .5 fim in order to draw a better comparison to 
the results of BO05. The largest, 2.5 yum sized particles show 
the broad plateau between 9/j.m and 12/j.m (see Fig. 3, two up- 
per panels) which does not appear to be existent in any of the 
observations modeled here. The difference in profiles between 
1.5/im and 2.5 /im sized particles from forsterite, enstatite and 
quartz is not significant (Fig. 3, three lower panels). Therefore, 
the 2.5 /im grains are not considered in the fitting procedure. 
Dust grains species which do not show an emission feature in 
the 8-13/im range, such as amorphous, compact carbon grains 



also known as mass absorption coefficient; e.g., Li 2005 
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Fig. 2. x 2 °f the fits versus the pre-specified temperature T. The ob- 
tained temperature T for the single fit and the object name are indi- 
cated in each plot. 



(e.g., Jager et al. 1998a), are not considered. However, consid- 
ering non-homogeneous dust grains (see Sect. 5), which consist 
of a mixture of various materials, carbon may also modify the 
silicate emission (see Natta et al. 2001). 

The emission profiles of the remaining dust components are 
linearly independent from the opacities of the other dust com- 
ponents as BO05 showed. This linear independence and the re- 
sulting uniqueness of the solution is not ensured if too many 
emission profiles are included. 

Previous investigations (e.g., Natta et al. 2000) do not in- 
clude certain dust species like enstatite, quartz and polycyclic 
aromatic hydrocarbons (PAH) in circumstellar disks around 
lower-mass stars. In general, this assumption seems not to be 



justified. Honda et al. (2003) found large mass fractions of en- 
statite and quartz in the TTauri object Hen 3-600 A. Finally, 
the robustness of our fitting procedure allows to include these 
dust species with one exception: we do not consider PAHs. The 
PAH features are strongly variable in peak position and profile 
from source to source and, therefore, can not be represented by 
one unique emission profile. Additionally, there are hints that 
the PAH features seem to be composed of several subfeatures 
(Peeters et al. 2002). 

4. Results 

In the following, we focus on the fitting results (see Table 2) 
and draw comparisons to previous publications. 

The results of the fitting procedure are satisfactory (see 
Fig. 4), in general. The quality of each fit, characterised by x 2 
(e.g., Press et al. 1986) with 



X 



1 N ^ 
J i=i 



N A -f 



cr v model (ii) - 3T m w) 
(TWw 



(3) 



is in the order of unity for most objects. The data with the 
largest^ 2 have the worst signal-to-noise ratios (e.g., SZ82, 
VZCha). In Eq. 3, Na is the number of data points, / is the de- 
gree of freedom, and Ti^dinot me measurement error of each 
spectral data point. Following PR03 and KS05 the measure- 
ment error is the standard deviation of spectra that was ob- 
served in subsequent nights. A root mean square (rms) of ~9% 
and ~11% were determined for TIMMI2 and LWS observa- 
tions, respectively. However, these errors were derived consid- 
ering also parts of the spectra with a significantly lower signal- 
to-noise ratio such as the boundaries of the N band and the 
spectral range where the telluric ozone-band appears. In the 
following, we assume a rms of ~6% and ~7%, respectively, 
which we derive from the wavelength interval between 10 fim 
and 12/iin. For several objects with higher signal-to-noise ra- 
tios, we determine a rms by using 



1 



N, 



NA-f 



J](rr del (A) 



■ r^\A)f. 



(4) 



i=i 



The errors on the fit parameters are calculated using a 
Monte Carlo method. For every spectrum we generate 100 syn- 
thetic spectra, by randomly adding Gaussian noise to the spec- 
trum at each wavelength point. The Gaussian noise has a rms 
of ~6% and ~7%, respectively. On each of these spectra, we 
perform the exact same compositional fit procedure, yielding 
different values of the fit parameters. From the resulting distri- 
bution of all fit parameters, we calculate the arithmetic mean 
value and standard deviation. 

Approximately at 9.7 /im a narrow "feature" appears in 
some spectra (e.g., Haro 1-16) which can not be fitted by the 
available emission profiles. We ascribe this artefact to a rem- 
nant of data reduction, i.e. sky subtraction, considering the 
spectral-adjacent, broad atmospheric ozone line at ~ 9.5 /xm. 
Additionally, there is a defect readout channel of the TIMMI2 
detector between 9.14/im and 9.67 /im. Because of these two 
issues, the fitting procedure does not account for the wave- 
length interval between ~9.2/mi and ~9.9^m, which is marked 
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by a dotted line in the Fig. 4. In some spectra (e.g., DQTau) a 
second remnant of sky subtraction appears at 12.7 //m, where 
several telluric OH-lines accumulate. In order to improve the fit 
quality, we additionally cut the spectral data in a narrow wave- 
length interval at ~ 12.7//m. 

Although PAH molecules appear to be abundant through- 
out the universe (e.g., Lagache et al. 2004; Wu et al. 2005), 
our spectra can be sufficiently fitted without considering PAHs 
although they also show an emission feature near 1 1 .2 yum si- 
milar to forsterite. The absence of its 11.3/im feature in the 
circumstellar disks of (our) T Tauri stars can be explained by 
the weak stellar UV radiation field compared to the higher- 
mass HAeBe stars. In fact, PAH molecules are thought to be 
excited by ultraviolet photons and cooled down by subsequent 
infrared emission. Furthermore, the circumstellar disk size and 
disk geometry affect the presence and strength of the PAH fea- 
ture as Habart et al. (2004) and Acke & van den Ancker (2004) 
showed. 

Pyroxene is an important part of our dust mixture because 
it is the dominant component at low temperatures (Gail 2003). 
BO05 included this dust component considering its contribu- 
tion to the short wavelength range of the silicate feature at 
~ 9/im (see Fig. 3). According to our results (see Table 2) 
the large mass fraction of large pyroxene grains (for example, 
in VZ Cha) is suspicious. A corresponding richness of large 
pyroxene grains was found in HAeBe stars (BO05), too. In 
1994, Pollack et al. determined that the initial composition of 
the silicate dust mixture in circumstellar disks only consists 
of a maximum of 30% pyroxene dust. Therefore, if the con- 
tribution of large pyroxene grains is real, small pyroxene dust 
grains have already disappeared due to grain growth. If we do 
not consider pyroxene 6 in our fitting function, the quantity x 2 
increases slightly, about a few percents for the T Tauri stars but 
up to 100% for the HAeBe stars HD 163296, HD 179218, and 
WW Vul. 

Apart from the objects Hen 3-600 A and S CrAN for most 
of the objects crystalline silicates contribute less than 10% to 
the total fit. This result agrees with Pollack et al. (1994) who 
assumed a crystalline initial contribution of only 10% while 
notable amounts of crystalline dust species like enstatite and 
forsterite are produced in annealing reactions during disk evo- 
lution. Gail (2003) shows that crystalline dust dominates in the 
inner parts of protoplanetary disks. However, the temperature is 
too high in the inner regions of the disk in order to produce the 
prominent silicate feature (see Fig. 1). Apparently, the crystalli- 
sation of dust is inhibited in the cool circumstellar disks around 
T Tauri stars which is confirmed by the minor mass fraction of 
crystalline dust in contrast to HAeBe stars where a mean crys- 
talline contribution of ~ 15.0% was determined (see BO05). 



Fig. 3. Opacities of grains with radii 0.1 yum (solid lines), 1.5 fim 
(dashed lines) and 2.5 fim (dotted lines). With increasing grain size 
the maxima shift to longer wavelengths and the emission profile of 
olivine and pyroxene flattens. 



large and small grains 



Table 2. Best results of our compositional fitting procedure. 



„ „, . t\k~\ Olivine Pyroxene Forsterite Enstatite Quartz 2 

Je ° L J Small [%] Large [%] Small [%] Large [%] Small [%] Large [%] Small [%] Large [%] Small [%] Large [%] X 



1 


GGTau 


384 ±3 


52.8 ± 2.3 


3.7 ±4.9 


1.2 ± 1.8 


34.6 ±5.1 


3.3 ± 0.2 


3.7 ±0.8 


0.2 ± 0.3 


0.4 ± 0.6 




- 


- 


0.88 ± 0.05 




2 


AATau 


442 ±6 


44.3 ± 5.6 


9.4 ± 11.8 


15.9 ±5.5 


23.9 ± 12.8 


- 


5.6 ± 1.1 


0.1 ±0.3 


0.2 ± 0.5 


0.1 


±0.2 


0.3 ± 0.5 


0.49 ± 0.02 




3 


LkCa 15 


1008 ± 95 


52.6 ± 3.3 


43.0 ±5.0 


0.1 ±0.5 


2.3 ± 3.3 


- 


0.4 ± 0.4 


1.0 ±0.4 


0.4 ± 0.6 




- 


0.2 ± 0.3 


0.70 ± 0.03 




4 


DQTau 


300 ±3 


53.3 ± 15.8 


3.1 ±5.3 


19.9 ± 13.2 


4.2 ± 6.8 


0.5 ± 0.4 


1.5 ± 1.2 


1.1 ± 1.2 


11.2 ±2.3 


2.0 


± 1.3 


3.2 ± 2.7 


0.62 ± 0.03 




5 


DRTau 


354 ± 1 


46.7 ± 3.3 


47.9 ± 4.3 


0.3 ± 0.9 


1.5 ±2.9 


1.1 ±0.3 


2.0 ± 0.8 


- 


0.2 ± 0.5 




- 


0.2 ± 0.4 


0.23 ± 0.01 




6 


GMAur 


403 ±6 


68.0 + 5.2 


27.5 ± 6.2 


- 


0.9 ±2.1 


1.2 ±0.3 


0.5 ± 0.6 


0.2 ± 0.3 


1.6 ± 1.1 




- 


0.1 ±0.3 


0.74 ± 0.03 


CO 


7 


SUAur 


348 ±2 


62.0+ 1.8 


18.9 ±4.1 


0.5 ± 1.0 


12.6 ±3.8 


2.9 ±0.1 


2.7 ± 0.5 


- 


0.2 ± 0.4 




- 


- 


0.55 ± 0.03 




8 


GWOri 


334 ±2 


54.0+ 1.1 


21.9 ±2.9 


0.4 ± 0.8 


22.1 ±2.8 


1.0 ±0.1 


0.3 ± 0.3 


- 


0.1 ±0.2 




- 


- 


0.37 ± 0.02 


'l'-i 
2 


9 


FUOri 


387 ± 1 


8.1 ±5.7 


82.7 ± 7.3 


0.2 ± 0.8 


2.4 ± 4.5 


3.5 ±0.6 


1.5 ± 1.5 


0.3 ± 0.6 


0.7 ± 1.3 




- 


0.6 ± 1.0 


0.20 ± 0.01 


a 


10 


BBW76 


324 ±3 


79.1 ± 8.6 


13.3 ± 8.0 


0.6 ± 1.7 


1.5 ±3.2 


0.2 ± 0.2 


1.1 ± 1.2 


0.2 ± 0.4 


0.5 ± 0.9 




- 


3.5 ± 1.7 


1.51 ±0.09 


EL 


11 


CRCha 


297 ±3 


59.7 ± 3.2 


2.9 ± 4.9 


7.0 ±2.9 


23.9 ± 5.7 


- 


5.0 ± 0.7 


0.2 ± 0.3 


1.0 ±0.8 




- 


0.2 ± 0.4 


2.31 ±0.42 




12 


TWHya 


285 ±3 


37.0 ± 4.4 


47.2 ± 7.4 


0.6 ± 1.5 


3.5 ±5.1 


2.9 ± 0.3 


5.9 ± 1.1 


- 


2.6 ± 1.5 




- 


0.2 ± 0.4 


3.32 ±0.58 


o 


13 


VWCha 


332 ±5 


3.1 ±4.2 


80.9 ± 11.3 


2.1 ±4.2 


8.1 ± 11.7 


0.2 ± 0.3 


0.7 ± 1.0 


0.5 ± 0.8 


0.5 ± 0.9 


0.6 


±0.4 


3.3 ± 1.4 


2.24 ±0.15 


3 


14 


Glass I 


453 + 2 


37.0 ± 0.7 


37.0 ± 2.0 


0.2 ± 0.4 


21.1 ± 1.9 


0.6 ±0.1 


1.4 ±0.2 


0.6 ± 0.2 


2.1 ±0.4 




- 


- 


0.05 ± 0.00 




15 


VZCha 


348 ±5 


4.0 ± 8.0 


13.8 ±21.9 


7.4 ± 11.7 


53.2 ± 24.3 


6.2 ± 1.7 


9.2 ± 4.6 


1.5 ±2.1 


1.6 ±2.7 


2.7 


± 1.0 


0.4 ± 0.8 


5.51 ±0.49 


C. 


16 


WWCha 


355 ±2 


14.2 ± 1.6 


35.9 ± 5.2 


1.2 ±2.0 


45.4 ± 5.6 


0.4 ±0.1 


1.4 ±0.4 


- 


0.1 ±0.2 




- 


1.4 ±0.3 


0.23 ± 0.01 


c 

:/■. 


17 


Hen3-600A 


301 ±4 


2.3 ± 4.7 


15.8 ± 19.6 


1.0 ±2.6 


36.6 ±21.2 


7.5 ± 1.2 


1.0 ± 1.5 


7.3 ± 3.0 


27.8 ± 3.7 


0.2 


±0.3 


0.4 ± 0.9 


0.70 ± 0.03 


-§. 


18 


IRAS 14050 


406 ±6 


74.9 ± 2.5 


2.4 ± 3.6 


0.5 ± 1.0 


19.1 ±4.1 


0.3 ± 0.2 


0.3 ± 0.4 


1.8 ±0.5 


0.6 ± 0.7 




- 


0.2 ± 0.4 


0.70 ± 0.03 


r. 


19 


SZ82 


411 ±7 


2.3 ± 3.6 


16.5 ± 15.8 


2.2 ± 4.4 


58.9 ± 16.8 


0.4 ± 0.4 


12.2 ± 2.4 


0.5 ± 0.9 


6.8 ± 2.5 


0.1 


±0.2 


0.2 ± 0.6 


1.52 ± 0.08 


c 

-i 
< 


20 


RULup 


346 ±2 


1.1 ±2.1 


5.5 ± 8.2 


2.7 ± 3.2 


75.0 ± 10.2 


6.9 ± 0.4 


0.3 ± 0.6 


3.5 ± 0.9 


1.8 ± 1.4 


0.7 


±0.4 


2.5 ± 1.1 


0.14 ±0.01 




21 


AS 205 N 


329 ± 1 


38.9 ±3.3 


5.7 ± 6.7 


2.0 ± 2.7 


47.5 ± 8.6 


1.9 ±0.3 


3.3 ± 0.9 


0.2 ± 0.3 


0.5 ± 0.8 




- 


- 


0.13 ±0.01 


c 


22 


AS 205 S 


507 ±4 


20.4 ± 1.3 


68.1 ±3.2 


0.5 ± 0.9 


9.5 ±3.1 


0.4 ±0.1 


0.9 ± 0.4 


- 


0.2 ± 0.3 




- 


- 


0.15 ±0.01 


H 

— ] 


23 


HBC639 


366 ±3 


2.3 ± 3.0 


69.8 ± 8.3 


0.5 ± 1.3 


13.5 ±9.3 


1.9 ± 0.4 


0.4 ± 0.5 


5.9 ± 1.0 


4.1 ± 1.7 




- 


1.4 ±0.7 


0.37 ± 0.01 


c 


24 


Harol-16 


373 ±4 


73.4 ± 2.5 


2.7 ±3.8 


1.0 ± 1.6 


18.9 ±5.9 


1.7 ± 0.3 


0.2 ± 0.4 


0.1 ±0.2 


1.8 ± 1.1 




- 


- 


0.69 ± 0.03 


V:. 


25 


AKSco 


559 ±5 


63.8 ± 1.1 


31.6 ± 1.7 


- 


0.5 ± 0.9 


1.7 ±0.1 


2.1 ±0.4 


- 


0.2 ± 0.3 




- 


- 


0.58 ± 0.03 




26 


SCrAN 


427 ± 1 


2.5 ±4.1 


42.2 ± 7.9 


- 


1.8 ±4.1 


18.4 ± 1.0 


11.6 ±2.3 


0.4 ± 0.7 


22.9 ± 2.6 




- 


0.3 ± 0.8 


0.15 ±0.01 


1 


27 


SCrAS 


375 ± 1 


1.0 ± 1.8 


2.7 ± 3.9 


1.0 ± 1.8 


89.1 ± 5.4 


3.8 ±0.3 


0.3 ± 0.5 


0.6 ± 0.6 


1.4 ± 1.1 




- 


0.1 ±0.2 


0.25 ± 0.01 




28 


MWC480 


801 ± 44 


0.7 ± 1.3 


59.2 ± 7.4 


0.1 ± 0.5 


33.9 ± 6.8 


1.2 ± 0.3 


3.0 ±0.9 


1.2 ±0.5 


0.5 ± 0.7 




- 


- 


0.69 ± 0.03 




29 


HD 163296 


431 ± 11 


20.2 ± 5.6 


24.7 ± 8.9 


0.8 ± 1.9 


44.5 ±5.1 


2.5 ± 0.3 


0.2 ± 0.4 


0.2 ± 0.3 


5.3 ± 1.2 




- 


1.6 ±0.6 


0.69 ± 0.02 




30 


HD 179218 


255 ±8 


1.2 ±2.1 


4.1 ±7.1 


1.7 ±3.7 


83.9 ± 7.6 


1.4 ± 0.2 


0.9 ± 0.7 


2.5 ± 1.0 


0.8 ± 0.9 


3.3 


±0.3 


0.2 ± 0.3 


0.69 ± 0.02 




31 


WWVul 


399 ±9 


52.7 ±5.1 


6.9 ± 7.2 


0.6 ± 1.3 


33.8 ±3.1 


1.6 ±0.2 


0.4 ± 0.5 


0.2 ± 0.3 


3.5 ± 1.0 




- 


0.3 ± 0.5 


0.69 ± 0.02 




32 


HD 184761 


519 + 10 


39.6 ± 3.4 


57.4 ± 4.0 


- 


0.4 ± 0.9 


0.7 ± 0.2 


1.6 ±0.7 


— 


0.1 ±0.3 




- 


- 


0.69 ± 0.02 





Note. The first two columns correspond to number (see Table 3) and name of each object. The value T represents the temperature of the underlying continuum which is determined assuming a black body function (see Eq. (2)). The abundance of 
small (0. 1 fim) and large (1.5 fim) grains of the various dust species are given as relative fractions of the total mass of silicate dust. We use the symbol '-' if a fitting parameter was not found in the corresponding spectra. The value x 2 refers to the 
quality of each fit (see Eq. (3)). 
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Analysing the results in Table 2, forsterite is the most abun- 
dant crystalline dust component in the circumstellar matter of 
TTauri systems in 19 of 32 objects. According to Gail (1998) 
forsterite is the thermodynamically most stable silicate com- 
pound and therefore the final product of the annealing proce- 
dure: amorphous grains — > polycrystalline grains — > enstatite 
— > forsterite. We point out that a back-reaction to enstatite 
is also possible if quartz is present: Mg2SiC»4 + SiCh — > 
2MgSiC>3 (Rietmeijer et al. 1986). Furthermore, the last step 
in this cascade can also be inhibited, e.g. by an O- or Mg- 
deficiency (Fabian et al. 2000) or by too low temperatures (Gail 
2003). Thus, the presence of enstatite in circumstellar disks of 
T Tauri stars can not be excluded although it appears less fre- 
quently than in HAeBe stars (BO05). 7 

Quartz is less abundant (see Table 2). This result is in agree- 
ment with the conclusion about the composition of the silicate 
mixture in stationary protoplanetary disks (see Table 8 in Gail 
2003). Gail (2003) determined a fraction of quartz which is 
< 5%. Furthermore, the small abundance of quartz -in contrast 
to HAeBe stars- could be a result of the quartz emission profile 
that we implement in our fitting procedure: our profile shows a 
more discontinous, less smooth curve in contrast to the emis- 
sion profile which BO05 used. 

We conclude by comparing our analysis with previous in- 
vestigations. Honda et al. (2003) succeeded to fit the mid- 
infrared spectra of the lOMyr old TTauri star Hen 3-600 A 
with larger amounts of quartz and less enstatite. The absence of 
quartz in our fitting results could be a result of "observational 
selection": the lower resolution of LWS (R=150) in compari- 
son with R=250 at the 8.2m Subaru Telescope equipped by the 
Cooled Mid-Infrared Camera and Spectrometer (Honda et al. 
2003). But the distinction can be also the consequence of the 
different fitting functions and the implementation of pyroxene 
and large crystalline grains in our fitting procedure. It is suspic- 
uous that our spectra tend to rise to longer wavelengths which 
indicates the negligible contribution of small amorphous grains 
in contrast to data of Honda et al. The modified shape of Hen 3- 
600 A could have its origin in a variable N band emission. 

N band spectra of the HAeBe stars HD 179218 and 
HD 163296 have been previously analysed by BO05, where 
a higher fraction of large enstatite grains was found for the 
object HD 179218. This distinction may be due to different 
spectral resolutions of the instruments (LWS versus TIMMI2) 
and differences in the fitting procedure. Our fitting results for 
HD 163296 confirm the results of BO05. 

4.1. Correlations 

While a correlation between feature strength and feature shape 
has already been observed in HAeBe and T Tauri stars indicat- 
ing grain growth (PR03; Acke & van den Ancker 2004; BO05), 
we examine the question if there is any correlation between 



grain growth, crystallinity and stellar properties like stellar age 
and luminosity. Here, we include only TTauri systems. 

1 . The degree of crystallinity in T Tauri systems correlates 
with neither stellar age nor luminosity (Fig. 7, 8). It is not 
clear if the lack of such correlations is a consequence of the 
uncertainty which adheres to the measurements in Table 
1. However, we suppose that a correlation 'crystallinity - 
stellar age' and 'crystallinity - stellar luminosity', respec- 
tively, exists in the inner circumstellar regions, where dust 
is more exposed to stellar irradiation. Interferometric obser- 
vations with the Mid-Infrared Interferometer (MIDI) show 
that the degree of crystallinity increases with decreasing ra- 
dial distances to the star (van Boekel et al. 2004). But the 
dust emission from these inner and hotter regions of the 
circumstellar disk (>700K) contributes only marginally to 
our spectra (Fig. 4 - 6) as it is illustrated in Fig. 1. 
With respect of HAeBe stars, BO05 did not investigate the 
correlations between crystallinity versus the stellar age and 
the stellar luminosity, respectively. Nevertheless, their ana- 
lysis of the correlation between the crystalline mass frac- 
tion and the stellar mass yielded a weak dependence: ob- 
jects with larger stellar masses tend to have larger mass 
fractions of crystalline grains. Assuming a correlation be- 
tween stellar mass and stellar luminosity this result dis- 
agrees with our study at first glance. But the correlation 
that was found by BO05 appears only for the more massive 
HAeBe stars (^ 3.0M o ). Those HAeBe stars which are less 
massive are more comparable to TTauri stars. 




7 We want to mention that recent mineralogical examinations have 
shown that enstatite, which could be developed from forsterite, is a 
final stage of silicate evolution in interplanetary medium and mete- 
orites, in particular, (see Fig. 12; private communication with Dominik 
Hezel, Institute for Mineralogy, Cologne). 



5 10 

stellar age [Myr] 



Fig. 7. Mass fraction in crystalline grains vs. stellar age according to 
Table 2. Numbers of objects are indicated. 

2. A clear correlation between the mass fraction of large 
grains (including large crystalline grains) and stellar lumi- 
nosity does not exist in our sample of TTauri stars (Fig. 
9). If we exclude objects with the worst signal-to-noise ra- 
tio Of 2 , TW Hya, VZ Cha and Sz 82) we find that most data 
points in Fig. 9 are located in the area which is limited by 
the dashed lines. The resulting, potentially existent correla- 
tion is weak because of the lack of objects with larger and 
smaller stellar luminosities (> 4L and < 0.5L©). BO05 
analysed the correlation between stellar mass and mass 
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stellor luminosity [!-,„„] 

Fig. 8. Mass fraction in crystalline grains vs. stellar luminosity. 
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Fig. 10. Mass fraction in large grains vs. stellar age. 
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Fig. 9. Mass fraction in large grains versus stellar luminosity. The stel- 
lar luminosity is plotted in logarithmic scale. Excluding objects with 
the worst signal-to-noise ratio most objects are located in the area 
which is limited by the dashed lines. The number of each data point 
corresponds to the object. 
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Fig. 11. Relative mass fraction in crystalline grains plotted versus the 
relative mass fraction in large grains. The number of each data point 
corresponds to the object in Table 2. 



fraction in large grains. Again, a correlation between both 
parameters exists only for more massive objects: HAeBe 
stars with a mass M k, 2.5M have a mass fraction of large 
grains above 85%. Less massive HAeBe stars do not show 
any dependence on the mass of the central star. 
We have to point out, that the stellar age does not affect the 
existence of large grains, which a further study shows (see 
Fig. 10). 

The Fig. 1 1 shows the relative mass fraction in large crys- 
talline grains versus the relative mass fraction in large 
grains. Apart from the objects Hen 3-600 A, SCrAN and 
DQ Tau most objects show a low mass fraction of crys- 
talline grains (<20%). Such a correlation could also be de- 
rived in HAeBe stars (BO05) where crystalline dust was ex- 
clusively found in objects with bigger amounts (>80%) of 
large grains. Around half of our T Tauri objects have only a 
negligible degree of crystallisation in spite of a high mass 
fraction in large grains. 



A correlation between grain growth of amorphous grains 
and grain growth of crystalline grains can not be found. 8 
The existence of such a correlation strongly depends on the 
rate of crystallisation and grain growth, respectively, and 
the ratio between both quantities. Note also that in some 
conditions the amorphisation of crystalline silicate occurs 
after ion irradiation (see, e.g., Djouadi et al. 2005). 
An agreement between T Tauri and HAeBe systems is ap- 
parent if the enstatite fraction of crystalline grains is plot- 
ted versus the mass fraction of crystalline grains. BO05 
found increasing enstatite fractions with increasing degrees 
of crystallinity. A similar tendency can be found in our sam- 
ple (see Fig. 12). The explanation for this weak correlation 
is a delayed production of forsterite or the reproduction of 
enstatite from forsterite (see Sect. 4, BO05 and Gail 2003). 



8 Both mass fractions are normalised by the mass fraction of crys- 
talline and the mass fraction of amorphous grains, respectively, con- 
sidering small and large grains, simultaneously. 
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Fig. 12. Enstatite fraction, normalised by the mass fraction of crys- 
talline grains plotted vs. the mass fraction of crystalline grains. 



Summarising our study, the analysis of the silicate feature of 
T Tauri systems is complicated by their generally lower bright- 
ness in contrast to HAeBe stars (BO05). However, we find cor- 
relations which partly agree with the analysis of BO05, for in- 
stance if their correlations are extrapolated to less massive sys- 
tems. 



5. Composite grains 

The opacities Kabs;i(/0 are conventionally calculated by Mie 
theory assuming homogeneous dust particles. In reality, col- 
lisions with gas, other grains and the interaction with radiation 
can lead to grain coagulation, partial grain evaporation, crys- 
tallisation, annealing, amorphisation, etc. As a result, compo- 
site, partially flaky and inhomogeneous dust grains with in- 
clusions of vacuum evolve (Dominik et al. 1993; Poppe et al. 
2000a,b). The calculations of the optical properties of such 
grains is a rather difficult problem which can be solved us- 
ing special algorithms like the discrete dipole approximation 
(DDA; e.g., Andersen et al. 2002, 2003). 

The most often used approximation includes the con- 
sideration of homogeneous particles with an average (effective) 
refractive index derived from effective medium theory (EMT; 
see Bruggeman 1935; Chylek et al. 2000). Comparing the DDA 
and EMT-Mie calculations, Voshchinnikov et al. (2005a) have 
recently discovered that the EMT-Mie approach can give rela- 
tively accurate results only for porous particles with "Rayleigh" 
inclusions. 9 Otherwise, the approach becomes unacceptable 
when the porosity exceeds ~50% of the grain volume. At the 
same time, the optical properties of heterogeneous spherical 
particles having inclusions of various sizes (Rayleigh and non- 
Rayleigh) and very large porosity were found to resemble well 
those of spheres with a large number (<; 15 - 20) of concentric 
layers with different compositions. 



It was also found that the layered-sphere model predicts a 
broadening of silicate bands and a shift of the peak position 
to larger wavelengths as porosity grows while in the case of 
the EMT-Mie model variations of the feature profile are less 
significant (see Voshchinnikov et al. 2005b for details). 

Since the cosmic particles can have inclusions of different 
sizes, here, we apply the model of multi-layered spheres for ex- 
planation of the 10/im silicate feature. For simplicity, we con- 
sider particles consisting of two materials (carbon and silicate) 
and 18 concentric layers 10 . Vacuum can be one of the com- 
ponents, too. In contrast to previous modelling with compact 
grains, the model parameters are grain size, composition and 
porosity. 

In the context of the described layered sphere model, pre- 
liminary calculations showed that variations of the volume 
fractions of vacuum, silicate and carbon produce diverse pro- 
files of the 10 /im feature from strongly peaked to very flat 
(Voshchinnikov et al. 2005b). It is important to note that differ- 
ent profiles are obtained for small grains with the same compact 
grain radius r cornp act = 0.1 yum. 11 This means that the variations 
of the shape of the silicate band can be related to the change of 
grain porosity and composition but not necessarily to the grain 
growth. 

As illustration Fig. 13 shows the result of our modelling for 
the objects AKSco and SCrAS having very different shapes 
of silicate band in their spectra. We used glassy olivine, va- 
cuum and amorphous carbon Bel (Rouleau & Martin 1991) 
and grains with a compact radius /" compact = 0.1 /im in both 
cases. Crystalline components are not considered. Figure 13 
and Table 3 show our best results after varying the relative 



Table 3. Parameters of composite grains used for modelling silicate 
band in spectra AK Sco and S CrA S (Fig. 13). 



object 



olivine 



Bel 



X 



AKSco 
SCrAS 



0.12yum 
0.27 yum 



49.5% 
3.35% 



5.5% 
1.65% 



45% 
95% 



0.57 
7.81 



Note. The quantity r por ous is the total radius of the porous grains (with vacuum). In col. 3-5 
volume fractions of materials are given. The parameter^ 2 is explained in Eq. 3. 



volume fractions of the three dust components. If the longer 
wavelength regime (A > 12.3 yum) of SCrAS is excluded the 
fits are satisfying according to Eq. (3). Considering this result 
and previous fitting with compact grains, the question about 
the unambiguous interpretation certainly arises. An answer on 
it can be found in the use of additional observations capable to 
discriminate between these two models (see Sect. 6). 

There are previous attempts to compare the silicate fea- 
ture with the mid-infrared emission behaviour of composite, 
complex dust grains. Bradley et al. (1999) found that the sub- 
micron components of interplanetary dust particles, glass with 
embedded metals and sulfides (GEMS), show features that are 
very similar to those in YSOs. The structure of GEMS and 



9 inclusions which are small in comparison with the wavelength of 
incident radiation 



10 Voshchinnikov et al. (2005a) have shown that this was enough to 
preclude influence of the order of materials on the results. 

1 ' The quantity r compact is the radius of the corresponding dust grain 
without voids. 
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Fig. 13. Normalised, continuum subtracted fluxes of AKSco and 
SCrAS. The dashed curves represent the results using the multi- 
layered model for porous grains. The relative volume fraction and the 
grain radius for AK Sco and S CrA S are presented in the Table 3. 

the properties of their iron-nickel inclusions refer to their ori- 
gin: Westphal et al. (2004) suggested that GEMS were prob- 
ably formed from crystalline grains by the exposure of hard, 
ionizing radiation in the stellar outflows from massive stars. 
Finally, the particles are accelerated in the shock waves of sub- 
sequent supernova explosions and distributed over space. The 
study of the importance of GEMS in circumstellar disks of 
YSOs and their actual influence on the 10/im-silicate-feature 
is out of the scope of this study. 

6. Conclusion 

In this study we investigated the mid-infrared silicate feature 
of 27 TTauri and 5 HAeBe stars. The relative mass fractions of 
amorphous and crystalline, large and small dust grains were de- 
rived by fitting a linear combination of the corresponding dust 
opacities. We assume that the emission feature has its origin in 
the optically thin surface layer of the circumstellar disk while 
optical thick disk layers only contribute to the underlying con- 
tinuum of the SED. Because of magnetorotational turbulences 
and turbulences which arise from Kelvin-Helmholtz instabil- 



ity dust grains with radii r < 1cm are assumed to be vertically 
well mixed in the disk (e.g., Johansen et al. 2005). The relative 
mass fractions obtained are a copy of the mean composition of 
micron-sized silicate grains in the whole disk, therefore. 

A central focus of this study was the question which phy- 
sical quantities can be derived from the observation profile of 
the silicate feature, in general. In terms of an impending vio- 
lation of the uniqueness of the fitting results the narrow spectral 
interval has a very restrictive effect on the number of fitting pa- 
rameters. Therefore, the number of the selected fitting param- 
eters should be as low as possible. Apart from the temperature 
T and a constant term Co we chose the opacities K a t, s of follow- 
ing dust species as fitting parameters: 0.1 fan- and 1.5/im-sized 
olivine, pyroxene, forsterite, enstatite and quartz. PAH emis- 
sivities were not implemented. The satisfying fits (see Table 2) 
underline the negligible abundance of excited PAHs in T Tauri 
systems in contrast to HAeBe systems. However, forsterite 
could overlay the PAH emission as it has a feature maxima at 
113fim, as well. 

As a subsequent result we found hints that both, grain 
growth and crystallisation appear simultaneously in (our sam- 
ple of) T Tauri systems while grain growth is more dominant. 
This result has to been seen in the context of previous investi- 
gations: Mathis et al. (1977) found a particle size distribution 
of several silicates in interstellar dust that can be reproduced by 
a power law of a strongly decreasing grain number with increa- 
sing grains size. In fact, Bouwman et al. (2001) found negligi- 
ble mass fractions of large silicate (> 1.5/im) grains in inter- 
stellar dust as parent material of dust in circumstellar (HAeBe) 
disks. However, the conclusion of grain growth is valid only 
for the case of compact dust grains which consist of one silicate 
component. Assuming inhomogeneous and flaky dust particles, 
which can be described by the multi-layered sphere model, for 
instance, the increase of porosity of composite grains up to 
50% and higher can be considered as an alternative to grain 
growth as it generates similar, flat spectra (Voshchinnikov et al. 
2005b). 

There are indications that all particles larger than 0.1 ^m 
are porous. But it is still a central question how large the de- 
gree of porosity is or if there is even a porosity distribution 
for grains of certain sizes. Furthermore, the EMT-model and 
the layered sphere model for porous grains are still competing 
models and partly provide different results. Although the lay- 
ered sphere model has its preferences (e.g., Rayleigh and non- 
Rayleigh inclusion can be handled) it focuses only on spherical 
grains. Grains with crystalline inclusion have not been consid- 
ered, so far. 

In respect of our fitting results, it is not clear, at all, how 
strongly porosity affects our fitting results. But according to 
the layered sphere model the contribution of large grains would 
be overestimated as a flat 10 //m feature can arise from large 
grains of any porosity or from very porous, small grains while 
a narrow feature results from small grains, exclusively. 

Considering the fact that T Tauri systems are less massive 
than HAeBe stars our analysis partly agrees with previous re- 
sults. In fact, our analysis showed that - in terms of the pro- 
perties of the circumstellar dust - T Tauri systems are a contin- 
uation of HAeBe systems at their lower mass end. Further ob- 
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servations - in particular, the silicate feature at 20 fim (Kessler- 
Silacci et al. 2006) - and polarization in silicate bands accom- 
panying with more comprehensive modelling could provide 
stronger constraints for existent silicate species in the circum- 
stellar disks of young stars. 

Certainly, the determination of the individual silicate com- 
position of the disk material is an important result, by its own, 
but it has also a further consequence for radiative transfer si- 
mulations. While recent radiative transfer simulations of cir- 
cumstellar disks used to be based on the extinction coefficients 
of an averaged mixture of silicates (e.g., "astronomical sili- 
cate", Draine & Lee 1984), the consideration of the results of 
compositional fits such as those presented here are meant to 
significantly improve disk modelling. Radiative transfer mod- 
els finally provide further constraints with respect of the fun- 
damental characteristics of the objects like geometry and struc- 
ture, sky position (inclination), (stellar) luminosity and object 
mass (star and disk). Together with the chemical composi- 
tion of the grains, models of YSOs become quite individual. 
Furthermore, such complete disk models open up new perspec- 
tives, e.g. the "compositional" analysis of absorption features 
which are strongly affected by disk structure and inclination. 
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Fig. 4. N band spectra of our sources (solid lines). The best fits to our spectra are plotted, too (dashed lines). We cut the wavelength interval 
~ 9. 1 fan to ~ 9.7 yum from the spectra (dotted line) taking into account the defect readout channel of the TIMMI2 detector and in order to 
exclude remnants of the telluric ozone corrections. 
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Fig. 5. Continuation of figure 4. 
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Fig. 6. Continuation of figure 4. The spectra of HAeBe stars are shown. 



